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Abstract 


Subscripts 


Results of an experimental investigation of a sym- 
metric crossing shock/turbulent boundary layer interac- 
tion are presented for a Mach number of 3.44 and deflec- 
tions angles of 2,6,8 and 9°. The interaction strengths 
vary from weak to strong enough to cause a large region 
of separated flow. Measured quantities include surface 
static pressure (both steady and unsteady) and flowfield 
Pitot pressures. Pitot profiles in the plane of symmetry 
through the interaction region are shown for various de- 
flection angles. Oil flow visualization and the results of a 
trace gas streamline tracking technique are also presented. 

Nomenclature 


So - 

D 

f 

G(f) - 

L.E. = 

Lgen = 

M 

P 

V 

Pt 

Pt2 = 

<7 = 

T.E. = 

v = 

x,y,z = 

x crs = 

x/e, zle = 

P 


shock generator plate deflection angle 

upstream B. L. thickness 

width and height of wind tunnel 

(30.5 cm) 

frequency (Hz) 

one-sided power spectrum 

leading edge 

length of shock generator plate 
Mach number 
static pressure 

fluctuating component of static pressure 
total pressure 
Pitot pressure 

RMS of fluctuating quantity 

trailing edge 

volume (mole) fraction 

cartesian coordinate system 

axial coordinate with origin at inviscid 

shock crossing location 

position coordinates of generator 

plate leading edge 

density 


'Aerospace Engineer, Inlet, Duct, and Nozzle Row 
Physics Branch, Member AIAA. 

Copyright © 1992 by the Am erica n Institute of Aeronautic* cod Astronautics, be. No 


w = wall condition 

0 = upstream reference condition 

oo = reference freestream condition 


Introduction 

O NE important class of shock/boundary-layer inter- 
actions is the glancing shock/turbulent boundary- 
layer interaction- This three-dimensional interaction oc- 
curs in many practical applications such as supersonic in- 
lets, nozzle flows and supersonic combustors. The strong 
pressure gradients present in the interaction region create 
complex three-dimensional flows and can lead to flow 
separation. In some cases, two glancing shocks cross 
and produce a local region of strong adverse pressure 
gradient. This investigation focuses on crossing shocks 
of equal strength (symmetric) interacting with a turbulent 
boundary layer. Shown schematically in Fig. 1, this ex- 
perimental configuration has the potential to provide good 
validation data for CFD methods. With a uniform in- 
coming boundary-layer, the boundary conditions are rel- 
atively simple and well defined. In addition, for equal 
shock strengths the flow has a plane of symmetry re- 
ducing the required calculation volume. By varying the 
shock strengths, the interaction can be varied from weak 
to strong enough to cause separation. 1 The separated 
cases provide a particularly interesting test for CFD val- 
idation since the separation is not fixed by the geometry. 



Fig. 1 Crossing shock wave/boundary-layer 
interaction with reference coordinates. 
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f As used in this paper, the term “separation” will refer 
to a distinct region of flow recirculation. 


Table 1 Summary of previous experimental work. 


Ref. 

Moo 

Re' /m 

Res 0 

6 0 mm 

a 

x crs /6 0 @T.E. 

Mee et al .' 

1.85 

3.9x 10 7 

7.4 xlO 4 

1.9 

5° 

3.8 

Batcho et 
al? 

2.95 

6.4 xlO 7 

2.4 x10 s 

4.0 

7° - 11° 

16-21 

Poddar et 
al? 

2.95 

6.6 xlO 7 

2.4 x10 s 

5.0 

7°, 11° 

16, 21 


2.47 

1.7 xlO 7 

65x10 s 

30.6 

4,6,8,9° 

2.4 - 2.8 

Hingst and 

2.96 

1.6x 10 7 

5.7 x10 s 

26.8 

4,6,8,10,11° 

1.6 - 2.5 

Williams 4 

3.44 

1.5 xlO 7 

5.1x10 s 

30.0 

4,6,8,10,12° 

0.6 - 1.7 


3.98 

1.3xl0 7 

4.7 x10 s 

33.6 

4,6,8,10,12° 

-0.3 - 1.0 

Garrison and 

2.96 

8.6 xlO 7 

2.6 x10 s 

3.0 

7° - 13° 

8.6 - 14.8 

Settles 6 

3.83 

7.9 xlO 7 

2.4 x10 s 

3.0 

o 

i 

o 

-0.2 - 10.5 

Kussoy and 
Horstman 7 

8.3 

5.3 xlO 7 

1.7x10 s 

32.5 

o 

wn 

o 

© 

1.0, 0.1 


Previous experimental investigations of the equal 
strength crossing shock/boundary-layer interaction are 
summarized in Thble 1. Some of these studies also con- 
sidered unequal strength crossing shocks. Mee et al? 
were the first to study the crossing shock interaction, 
providing mean surface static pressure data and limit- 
ing surface streamline patterns. Batcho et al 2 provided 
similar data but also included preliminary measurements 
of the unsteady surface pressure and temperature along 
the plane of symmetry. These results showed that, un- 
like the data of Mee, for all deflection angles considered, 
the pressure rise downstream of the inviscid shock cross- 
ing location significantly exceeded (by as much as 50%) 
the theoretical inviscid shock pressure rise. Using the 
same facility as Batcho, Poddar et al? mapped the un- 
steady surface pressure field. These measurements indi- 
cate that for stronger interactions, the flow downstream of 
the inviscid shock crossing position has a significant un- 
steady characteristic. Hingst and Williams 4 * 5 presented 
surface oil flow patterns and mean surface static pres- 
sure distributions over a relatively wide range of oper- 
ating conditions. The first flowfield data for the cross- 
ing shock configuration was obtained by Garrison and 
Settles 6 who performed laser sheet flow visualization in 
the cross plane. From these results, a model of the com- 
plex shock structure through the interaction region was 
constructed. Most recently, Kussoy and Hors tm an 7 per- 
formed surface flow visualization, surface heat transfer 
and Pitot pressure measurements in the cross plane for 
hypersonic flow conditions. 

Several of the experimental data sets have been 
used as a comparison to numerically generated results. 
Gaitonde and Knight 8 calculated the Mach 1.85, a= 5° in- 
teraction of Mee et al. Qualitatively, agreement with the 
experiment was good, but the calculations significantly 
overpredicted upstream influence and the rate of pres- 


sure rise along the plane of symmetry. Later, Gaitonde 
and Knight 9 numerically investigated a stronger symmet- 
ric crossing shock configuration (M=3, a =4° and 8°) 
as a means for controlling boundary layer separation, 
but no experimental data were available for compari- 
son. Narayanswami et al 10 performed calculations based 
on the experimental configuration of Batcho and Poddar. 
Calculations were performed for M=2.95, a=9° and 11° 
utilizing both the algebraic turbulence model of Baldwin 
and Lomax 11 and the k-e model of Jones and Launder 12 . 
Overall, the calculations with the algebraic model agreed 
fairly well with the experimental surface pressure data, 
but upstream influence was underpredicted. The k-e cal- 
culations were considered preliminary and required fur- 
ther refinement The data of Hingst and Williams 4 * 5 were 
recently used as a comparison for CFD calculations per- 
formed by Reddy 13 using the PARC3D code in conjunc- 
tion with the Baldwin and Lomax turbulence model. The 
M=3.5, a =6° and 10° and M=4.0, a =6° cases were com- 
puted. The results of these calculations showed overall 
good agreement with the surface pressure data for both 
the unseparated and separated cases. For the separated 
case, however, the reverse flow region was predicted to 
be further upstream and larger than that observed exper- 
imentally. 

A common feature of the crossing shock experiments 
is that, at the trailing edge of the generator plates, the 
flow encounters an expansion. The question then arises 
as to what degree does the upstream influence of the trail- 
ing edge expansion affect the interaction region. Clearly, 
if the expansion does feed upstream into the interaction 
region, it will tend to relieve the streamwise pressure 
gradient and thus the strength of the interaction. The dis- 
tance between the shock crossing location and the gen- 
erator plate trailing edge plane is a relevant parameter 
with regard to downstream influences. The last column 
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in Table 1 enumerates this distance normalized by the 
incoming boundary layer thickness for the various exper- 
imental studies. Of the studies where this parameter is 
relatively small, measured centerline static pressure was 
never observed to exceed that which is predicted by invis- 
cid theory. The relatively large value of this parameter 
for the configuration tested by Batcho and Poddar cor- 
relates with the large centerline pressure rise above the 
inviscid value that they observed. One of the goals of this 
study is to examine the effect of varying this parameter 
by changing the generator plate length. 

Experimental Program 


Test Configuration 

The present investigation was conducted in the 
NASA Lewis Research Center 1 x 1 ft Supersonic Wind 
Tunnel. This wind tunnel is a continuous flow facility 
with Mach number variation provided by interchangeable 
nozzle blocks. The crossing shock/boundary-layer exper- 
iment is configured by using two shock generator plates 
that span the wind tunnel test section. The shock gener- 
ators, when at angle of attack to the free stream tunnel 
flow, produce oblique shocks. The intersection of these 
shocks with the naturally occurring boundary-layer on the 
tunnel walls defines the experiment. A schematic of the 
test configuration with reference coordinates is shown in 
Fig. 1. With the exception of the shock generator plate 
lengths, the test hardware and actuation system for the 
present study is the same as was used by Hingst and 
Williams. 4 Deflection of the plates is accomplished by 
rotating the plates about a point For reference, the posi- 
tion of the leading edge relative to the upstream reference 
plane (x=0) and the wind tunnel centerline (z=0) is given 
by the following (see Fig. 1): 

i/e(cm) = 6.350 — (9.525sin a + 3.810 cosa) (1) 

:/e(cm) = 15.875 — (9.525 cos a — 3.810 sin a) (2) 

Two sets of generator plates which are longer (L gcn = 
22.9 cm and 25.4 cm) than the ones used by Hingst and 
Williams (L gcn = 20.3 cm) were fabricated in order to 
investigate how generator length affects the interaction 
region. The additional length is added to the generator 
plate trailing edge so that eqs. (1) and (2) are applicable 
for all plate lengths. 

Results and Discussion 

The initial plans for the test called for Moo =3 .44 and 
a=6° and 10°. These conditions were previously studied 
experimentally by Hingst and Williams 4 and numerically 
by Reddy 13 and correspond to unseparated and separated 
flow conditions, respectively. A significant difference, 
however, was that longer generator plates would be used 


in the present study to reduce the upstream influence 
of the trailing edge. To this end, the wind tunnel was 
initially configured with the 25.4 cm generator plates 
(the aforementioned studies considered 20.3 cm generator 
plates). At Mach 3.44, these plates were found to be good 
only up to deflection angles of 9°, after which the flow 
between the plates would completely unstarL With these 
plates, only surface static pressure and oil flow patterns 
were obtained for deflection angles between 2° and 9°. 
With the 22.9 cm plates installed, deflection angles up to 
11° could be attained without a probe and up to 9° with a 
probe. It was then decided that the 22.9 cm plates would 
be used for flow field measurements, recognizing that 9° 
is the upper limit for deflection angle. Some limited data 
were also obtained at Moo =2.96 for which the maximum 
deflection angle with a probe in the flow was determined 
to be 8° , but none is included in the present report 

Upstream Flow Conditions 

The flowfield upstream of the generator plates was 
surveyed with a Pitot probe to insure that wind tunnel 
comer effects and the distorted boundary layers, which 
develop on the non-contoured wind tunnel nozzle walls, 
were bypassed behind the generator plates. Fig. 2 shows 
the results of the surveys plotted in terms of Mach number 
contours. The left and right sides of this figure correspond 
to M=2.96 and 3.44, respectively. Also indicated is the 
position of the shock generators at the largest deflection 
angle considered for each freestream condition. The up- 
stream boundary layer was surveyed on the wind tunnel 
centerline (z=0) using a flattened Pitot probe. Various 
computed boundary layer parameters for the Mach 3.44 
case are summarized in Table 2. The skin friction coeffi- 
cient shown in the table was deduced from a curve fit of 
the data to the Sun-Childs 14 wall-wake profile. Analysis 
of the profile indicates that the boundary layer, which ef- 
fectively begins at the wind tunnel nozzle throat and is not 
subjected to cross flow pressure gradients in the nozzle, 
can be characterized as uniform and in equilibrium. 


Table 2 Upstream reference conditions (x=0). 


Moo 

3.44 

Pi, co 

241.3 

P w ,o kPa 

3.46 

T t , oo K 

297.0 

Re' /m 

1.5x 10 7 

Re 6o 

5.1x10 s 

So mm 

30.0 

mm 

10.7 

60 mm 

1.4 

Cj , 0 x 10 3 

1.18 
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Fig. 2 Mach number contours 
at upstream station (x=0). 

Flow Visualization 

Near-wall limiting streamline behavior was investi- 
gated by means of oil flow visualization. A powdered 


fluorescent dye is mixed with 140 wt. oil and painted 
on the surface between the generator plates. The wind 
tunnel is then run until a steady pattern is established and 
then rapidly shut down to preserve it Oil flow results 
were obtained for both the 22.9 cm and 25.4 cm plates. 
Three cases of strong interaction, i.e. the flow is sepa- 
rated, are shown Fig. 3. Fig. 3b shows the flow pattern 
for a deflection angle of 9° and a generator plate length 
of 22.9 cm. A significant separation region is clearly 
present Other oil flow tests, which are not presented 
here, suggest that incipient separation occurs at approxi- 
mately q=S° for the 22.9 cm long generator plates. Com- 
parison between Fig. 3b and Fig. 3a shows the effect of 
increasing the deflection angle from 9° to 10°. The sep- 
aration region increases in size and moves forward. An 
increase in deflection angle beyond 10° causes a com- 
plete unstart between the plates. Comparison between 
Fig. 3b and Fig. 3c illustrates the effect of increasing the 
plate length by 2.54 cm while maintaining the deflection 
angle at 9°. Although the strength of the shocks remain 
the same, the additional compression associated with the 
longer plates causes the separated region to move for- 
ward and grow in size. 



a) a=10°, Lgen = 22.9 cm. 

Fig. 3 Surface oil flow patterns. (Continued . . . ) 


4 



b) a=9°, Lo«„ = 22.9 cm. 



c) a =9°, Lgen = 25.4 cm. 

Fig. 3 Surface oil flow patterns. 
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Mean Surface Static Pressure 

Measurements of the mean surface static pressure 
over the entire surface between the generator plates were 
obtained using a repositionable plate with eighty 0.508 
mm diameter taps. The plate could be installed in four 
different positions which yields a total of 320 taps. How- 
ever, data from taps located under or behind the generator 
plates were discarded The location of the pressure taps 
relative to the generator plates are shown in Fig. 2 of 
Ref. 5. 

Contours of the surface static pressure for a Mach 
number of 3.44 and deflection angles between 2° and 10° 
are shown in Fig. 4. These data were obtained with the 
22.9 cm plates with the exception of the lower right plot 
which corresponds to the 25.4 cm plates. The shaded re- 
gions indicate the approximate region of flow separation 
as deduced from the oil flow tests. The “x” symbols 
on the centerline indicate the positions of Kulite dynamic 
pressure transducers. From the results shown in Fig. 4, it 
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INVISCID SHOCKS 


3 - X - KULITE LOCATION " 



a=8 deg., Lgen=22.9 cm - 



a=6 deg., Lgen=22.9 cm - 



-3 - o=2 deg., Lgen=22.9 cm " 


is apparent that there is a qualitative difference between 
the pressure distributions for the unseparated and sepa- 
rated cases. Considering the surface pressure distribution 
along transverse lines (x=constant), it is seen that in the 
initial part of the interaction region for both cases, the 
minimum pressure occurs on the centerline (z=0). At a 
point upstream of the shock crossing location, a transition 
occurs whereby the centerline pressure becomes a max- 
imum. For the unseparated cases, this behavior persists 
throughout the remainder of the interaction region. For 
the separated cases, however, a second transition occurs 
midway through the interaction such that the centerline 
pressure again becomes a minimum. This second transi- 
tion, which is characterized by a transverse line of nearly 
constant pressure between the generator plates, correlates 
well with the location of the separation region. An ex- 
amination of the oil flow and surface pressure distribu- 
tion data for both the unseparated and separated cases 
of Hingst and Williams 4 reveals the same behavior for 



a=10deg., Lgen=22.9cm - 



a=9 deg., Lgen=22.9 cm - 



-6 -4 -2 0 2 4 
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Fig. 4 Wall static pressure distributions. 
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the 20.3 cm plates. This distinguishing behavior between 
unseparated and separated flow conditions is also evident 
in the numerical calculations performed by Gaitonde and 
Knight 9 Also, a comparison of the surface flow visual- 
ization and surface pressure presented by Batcho el al. 2 , 
which for all cases were unseparated, agree qualitatively 
with the present unseparated cases. This observed differ- 
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ence in surface pressure behavior for the unseparated and 
separated cases may be useful as a means of detecting 
separation in symmetric crossing shock configurations. 

Fig. 5 shows the surface pressure distributions in 
the plane of symmetry for both the 22.9 cm and 25.4 
cm generator plates. Also, Hingst and Williams 4 data 
obtained with the 203 cm generator plates at Mach 3.44 



Xcry<5 0 


Fig. 5 Centerline wall static pressure distributions. 
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and a=4,6,8 and 10° are included for comparison. Note 
the scale change between a =4° and a=6 ° . At the weakest 
interaction, <*=2°, there is little difference between the 
distributions corresponding to the two longest generator 
plates. At a=4°, data for the two longer plates again 
agree well with each other. However, the distribution 
corresponding to the 20.3 cm plates agrees with the longer 
plate data only for a short distance into the interaction 
region, after which the short plate data lie below the 
longer plate data. At a=6 ° , a similar trend is observed for 
the shortest plate, but also the distributions for the longer 
plates begin to deviate. Further increases in the deflection 
angle results in larger deviations for all plate lengths. 
We may conclude from these results that even for small 
deflection angles, the interaction region is influenced by 
the trailing edge expansion when the 20.3 cm plates are 
used. Indeed, even the longer plates appear to capture the 
true strength of the interaction only for the weaker cases. 
Achieving total isolation from the trailing edge expansion 
in a given facility may be difficult To do so requires 
a relatively thin boundary layer which makes flow field 
measurements difficult The presence of a trailing edge 
expansion does not invalidate a data set, but does require 
that it be modelled accurately in a CFD calculation. 

Pitot Pressure Profiles in the Plane of Symmetry 

Axial development of Pitot pressure profiles in the 
plane of symmetry are shown for a= 6, 8 and 9° in 
Figs. 6, 7 and 8, respectively. For these measurements, 
a single flattened Pitot probe was used with outer tip di- 
mensions of 0.90 x 0.41 mm. The probe was inserted 
into the flow parallel to the generator plates from the wall 
opposing the interaction surface. Recall that all flowfield 
measurements were obtained with the 22.9 cm generator 
plates. The lower abscissa of the plots gives the axial lo- 
cation of the surveys relative to the inviscid shock cross- 
ing location and is normalized by the incoming boundary 
layer thickness. The upper abscissa is the survey location 
(in inches) measured from the upstream reference plane 
(see Fig. 1). At each axial profile station, the theoreti- 
cal inviscid Pitot pressure is shown for comparison with 
the measured values. Scaling for the Pitot pressure is 
indicated (in kPa) above the upperleftmost profile. For 
reference, the wall static pressure is indicated on the plots 
in two ways. First, it is plotted as a limiting value to the 
Pitot surveys as the wall is approached (solid symbols), 
and second, the axial distribution of wall static pressure 

is overlaid ( line) and scaled as indicated by the 

right hand ordinate. Also indicated in these plots are the 
locations of the Kuiite pressure transducers and the lo- 
cation of the trailing edge of the generator plates (T.E.). 
For a=9° (Fig. 8) the downstream region of the interac- 


tion could not be surveyed due to probe blockage caus- 
ing the model to unstart. It should also be noted that the 
presence of the probe may affect the flowfield when the 
flow is separated or at near separation conditions. These 
effects, if any, are difficult to quantify. 

Although the surveys were intended to pass directly 
from Zone 1 to Zone 3 through the shock crossing loca- 
tion (see Fig. 1), the results indicate that the actual axial 
traverse passed through a distance within Zone 2. This 
can be seen by examining the profiles immediately up- 
stream and downstream of the inviscid shock crossing 
location. In the core region of the profiles, the Pitot pres- 
sure does not agree with either the Zone 1 or the Zone 3 
theoretical values. In fact, they agree very well with the 
theoretical pressure in Zone 2. The axial extent of this 
deviation decreases with increasing deflection angle. The 
cause for the deviation is unknown, but may be related, 
at least in part, to probe deflection in the vicinity of the 
shock waves. 

With reference to Fig. 6 (a=6° configuration), the 
shock induced pressure gradient creates two distinct fea- 
tures in the profiles. First, the deceleration of the near- 
wall flow produces an inflection point in the profiles, and 
secondly, this low momentum near-wall fluid acts as an 
obstacle to the upstream flow, spawning an oblique shock 
wave. The presence of the shock is embodied in a local 
rise in Pitot pressure above the core value. The origin 
of the shock is a coalescing of compression waves in the 
boundary layer and appears as a smooth bump rather than 
a discrete discontinuity. As it rises out of the shear layer, 
it becomes a more discernible step (see e.g., the profile at 
x=7.0 in). In Garrison and Settles’ 6 model of the interac- 
tion, this shock is identified as the “centerline segment of 
the reflected separation shock” (segment “6a” in Fig. 12 
of Ref. 6). For all profiles, the transverse gradient of the 
Pilot pressure at the wall always remains greater than zero 
(dP t 2 /dy > 0) indicating that the flow does not separate. 
Increasing the deflection angle to a=8° (Fig. 7) produces 
similar results, but now the presence of the compression 
can be distinguished further upstream. Like the a=6° 
case, the gradient at the wall, although reduced, is non- 
zero throughout the interaction region. For a=9° (Fig. 8), 
the presence of the compression is observed even further 
upstream. The pressure gradient at the wall for this case 
appears to reach a zero value through a large portion of 
the interaction implying that the flow has separated. This 
region of zero gradient correlates well with the separated 
region observed in the oil flow results which are indi- 
cated in the figure by the shaded region. The step in the 
last axial survey is believed to be caused by an incipient 
unstart condition in the tunnel inasmuch as moving the 
probe further downstream resulted in a complete unstart 
of the model. 
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Fluctuating Surface Pressure 

Limited unsteady surface pressure data in the plane 
of symmetry were measured with Kulite Model XCS- 
062 dynamic pressure transducers. Four transducers were 
mounted on a repositionable plate that could be installed 
in four different positions. During data reduction, it was 
determined that the mean output of one transducer was 
erratic and hence not used. Thus, at a given operating 
condition, pressure traces were recorded at a total of 
twelve stream wise locations on the centerline. These 
locations are indicated by the x’s in Figs. 4, 6, 7 and 
8. Unfortunately, the removal of the faulty transducer 
resulted in a rather large axial gap in the interaction 
region. All the Kulite data were obtained with the 22.9 
cm generator plates. 

The output from a transducer was amplified by a 
Preston amplifier Model 8300 XWB with a gain setting 
of 100. The amplified signal was then lowpass filtered at 
50 kHz. The output of the filters was further amplified 
and offset via a Thermal Systems Inc. Model 157 signal 
conditioner. The gain and offset of the signal conditioner o 
were adjusted for each data record so as to take full 
advantage of the ±5 volt input range of the A/D board. 

The signal conditioner output was digitized and recorded * 
on a Concurrent Corporation data acquisition computer 
Model 5240. For all data records, the signal was sampled 
at 200 kHz/channel. The number of samples per channel 
was 204,800 resulting in a total sampling time of just over 
a second. A comparison between the mean Kulite output 
and the conventional wall pressure taps for a deflection 
angle of a=9° is shown in Fig. 9. 

Wall static pressure RMS distributions in the plane 
of symmetry are shown in Fig. 10a normalized by up- 
stream reference values and in Fig. 10b normalized by 
the local wall pressure. The left ordinate in in Fig. 10a 
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Fig. 9 Comparison between Kulite and wall 
tap mean pressure, a =9°, L gen =22.9 cm. 


is referenced to the undisturbed mean wall pressure and 
the right ordinate is referenced to the undisturbed RMS 
value. The upstream portion of the a=2° distributions are 
sufficiently removed from the interaction region so as to 
be representative of the undisturbed incoming boundary 
layer. The fluctuation intensity in this region is just under 
1%. For deflection angles greater than 2°, the distribu- 
tions exhibit a peak at the beginning of the interaction 
that coincides with the start of the mean pressure rise. 
For the strongest interaction considered (a=10°), the in- 
tensity rises very abruptly to a relatively high value. With 
reference to the upper right plot of Fig. 4, this high in- 
tensity was measured at the first Kulite location upstream 
of the separated region. The initial peak is followed by 
a decrease and then a second peak whose axial position 
varies with the deflection angle. 



a) Normalized by undisturbed quantities. 



X c J 6 0 

b) Normalized by local wall pressure. 

Fig. 10 Wall pressure RMS distributions 
in the plane of symmetry, Lge n =22.9 cm. 
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The results in Fig. 10 may be compared with the 
Mach 2.96 data of Batcho et al } Recall from Table 1 
that the boundary layer thickness for this study is nearly 
an order of magnitude thinner than the present case. A 
comparison of the RMS pressure levels for the undis- 
turbed boundary layer shows that they are virtually the 
same. Although the operating conditions are somewhat 
different, the a-T and a=ll° cases in Batcho’s study 
are qualitatively similar to the a=6°and o=9° cases, re- 
spectively, of the present study. Quantitatively, the RMS 
levels for Batcho’s a=T case fall between the present 
a= 6° and a=8° cases. 

Normalized power spectral density functions corre- 
sponding to each of the data points in Fig. 10 were cal- 
culated and are shown in Fig. 11. Along the top of the 
plots, the stream wise location of the Kulite transducer 
relative to the upstream reference plane is denoted. The 
number in parentheses is a reference number used to dis- 
tinguish between the three good transducers (transducer 
#3 was determined to be faulty). The transducers located 
at x=2.18, 2.81 and 3.31 inches for the a=2° case should 
be representative of the undisturbed incoming boundary 
layer. Each of these locations correspond to a different 
transducer. A notable difference between the three trans- 
ducer outputs is the presence of several high frequency 
(f > 10kHz) spikes for transducers #2 and #4 which are 
not present for transducer #1. The source of this high 
frequency noise has been traced to individual channels of 
the data acquisition A/D board. Another feature common 
to all three transducers in the undisturbed spectrum plots 
is the presence of a large amplitude fluctuation at approxi- 
mately 300 Hz. The source of this frequency is unknown, 
but some recent measurements made by Barnhart 15 in the 
NASA Lewis lxl SWT may shed some light as to their 
origin. As a component of his study of a blunt fin interac- 
tion, Barnhart recorded unsteady wall pressure traces for 
Mach 2, 3, 4 and 5 in the undisturbed boundary layer. 
Identical instrumentation was used with the exception of 
the second amplifier which in his study was a Pacific In- 
struments signal conditioner. For the Mach 2 case, Barn- 
hart’s data showed a very discrete spike at about 1kHz 
which, interestingly, all but disappeared when the tunnel 
mass flow rate was reduced. The second amplifier used 
in the present study was subsequently tested to determine 
if it was the source of the 300 Hz component It was 
not If this noise is truly a feature of the incoming flow 
it is puzzling as to what physically could create such a 
discrete frequency. Clearly, this issue must be resolved 
before any further measurements of this type are made. 


Ignoring for the moment the noise components of the 
signals, the frequency content of the incoming boundary 
layer is fairly broadband with the dominant frequency be- 
tween 2 to 3 kHz. With the exception of the a=10° case, 
the spectrum for the Kulite transducer located furthest up- 
stream (x=2.18 in) is independent of the generator plate 
deflection angle. At a=10°, some influence of the sep- 
aration shock is distinguishable in the spectrum. As the 
strength of the interaction is increased, the contribution of 
the noise components is diminished until at the strongest 
interaction considered (a=10°), they are negligible. For 
the a=6, 8, and 9° cases, moving through the interac- 
tion results in a slight shift in center frequency from the 
undisturbed value of ~2.5 kHz down to ~1.5 kHz. At 
the beginning of the or=10° interaction, there is a large 
shift in center frequency down to —150 Hz (see x=2.81 
in). This location corresponds to where the highest RMS 
fluctuations were observed and is the first Kulite location 
upstream from the separated region (see upper right plot 
in Fig. 4). The 150 Hz frequency can also be observed 
in the spectrum immediately downstream (x=331), but to 
a lessor degree. A similar low frequency component is 
also observed in the first data point upstream of the invis- 
cid crossing shock location for Poddar’s a=7° symmetric 
interaction. Surprisingly, this low frequency component 
was absent in Poddar’s a=ll° data, but this may be due 
to inadequate resolution of the region. Moving through 
the remainder of the interaction region, the spectrum for 
the a=10° case approaches a similar shape as the weaker 
interaction cases with one notable difference. The first 
pressure trace downstream of the inviscid shock cross- 
ing location (x=7.43 in) exhibits a second peak in the 
spectrum at about 10 kHz. A trace of this frequency 
component can also be seen in the spectrum immediately 
downstream (x=8.06 in). Unfortunately, the rather large 
axial gap in the data prevents complete assessment of 
this feature. Another feature of these results is the broad- 
ening of the spectrum as the shock crossing location is 
approached (see e.g., the distributions at x=4.93)., but 
again analysis is limited by the axial gap in the data. 

Through most of the interaction region of the present 
study, the dominant frequencies have been observed to be 
centered in the range from 1 to 3 kHz, excluding the low 
frequency peak observed at the beginning of the a=10° 
interaction. In comparison, the Mach 3, thin boundary 
layer results of Poddar et al? show dominant frequencies 
an order of magnitude higher. This difference is presum- 
ably related to the different length scales associated with 
the different boundary layer thicknesses. 
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Fig. 11 Wall pressure power spectral densities in the plane of symmetry. (Continued . 
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Fig. 11 Wall pressure power spectral densities in the plane of symmetry. 


Pitot Pressure Contours 

Pitot pressure surveys were conducted in the cross 
plane just upstream of the generator plate trailing edge 
plane. These measurements were made in a similar man- 
ner to the centerline profile surveys except that a flat- 
tened 3-hole Cobra probe was utilized. The measure- 
ment plane relative to the model geometry is shown in 
Fig. 12. Results of the surveys for a=6, 8 and 9° are 
shown in Fig. 13. In these plots, the local Pitot pressure 
is normalized by the undisturbed upstream value (P^o = 
55 kPa). Probing at this location for the a=9° case was 
possible because the main body of the probe was in the 
expansion region downstream of the trailing edge. These 
dam were taken in only one half of the facility and then 
imaged for presentation purposes. 


x - TRACE GAS INJECTION LOCATION 



Fig. 12 Transverse measurement plane. 



z/<5 0 


a) a = 6°. 

Fig. 13 Pitot pressure distributions (Ptt/Pt 2 ,o) at x 
= 9.25 in, Lgen = 22.9 cm. (Continued . . . ) 



- 2-10 1 2 
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b) a = 8°. 



z/<5o 

c) a = 9°. 


Fig. 13 Pitot pressure distributions 
(Pt 2 /Pt 2 ,o) at x = 9.25 in, Lge„ = 22.9 cm. 

In the core region of the flow, the measured shock 
location agrees fairly well with the inviscid theory. As 
the interaction region is approached from the core region 
(-y direction), the distance between the shocks increases. 
The lifting of the boundary layer leading to separation 
is clearly seen as the strength of the interaction is in- 
creased. The distortion of the contours is indicative of a 
pair of contra-rotating vortices centered about the plane 
of symmetry with common flow away from the floor sur- 
face. The presence of these vortices appears in previous 
calculations of the crossing shock configuration. 10 
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Volume Fraction Contours 

Near wall transport was investigated by means of 
a high speed trace gas technique recently developed by 
Reichert 16 and extended to supersonic flows by Davis 
and Hingst. 17 . The technique involves leaking pure ethy- 
lene through a static pressure port and then sampling the 
flow in a downstream plane. Sampling is done through 
a conventional Pitot tube and pumped into a flame ion- 
ization detector (Gow-Mac Model 23-500) hydrocarbon 
analyzer for analysis. The output from the analyzer is the 
local concentration of the ethylene-air mixture in ppm of 
ethylene. Ethylene was injected into the flow field from 
two locations. The location of the taps and the location 
of the measurement plane are shown in Fig. 12. Results 
for various deflection angles are shown for the centerline 
and offset injection in Fig. 14 and Fig. 15, respectively. 
For the centerline injection case, data were taken on only 
one side of the centerline and then imaged for presen- 
tation. The results are plotted in terms of the volume 
fraction normalized by the maximum volume fraction in 
the measurement plane for the a = 2° centerline injection 
case. The development of the ethylene plume is a func- 
tion of both mean convection and turbulent diffusion. For 
the or=2° centerline injection case, the peak concentration 
occurs on the surface and the remaining distribution is pri- 
marily the result of turbulent diffusion. At a =6° , the peak 



. M = 3.44 


a = 2.0 



Fig. 14 Volume fraction distributions for M = 3.5, 
a = 2°, 6° and 9° cases, centerline injection. 


concentration decreases by about 20% and has lifted from 
the surface by mean convection. The transverse width of 
the plume does not increase significantly over the a=2° 
case. At a =9°, the peak concentration decreases signif- 
icantly and is lifted even further from the surface. The 
width of the plume has increased considerably indicat- 
ing an substantial increase in turbulence intensity. An 
interesting aspect of these results is that the distributions 
do not have the characteristic horseshoe shape associated 
with the stream wise vortex structure that is implied by 
the Pitot pressure distributions. 

Similar comments also apply to the offset injection 
case (Fig. 15), but here a greater degree of distortion of 
the plume is observed, particularly for the a=9° where the 
elongation in the y direction implies an upward flow along 
the generator plates. Fluid for the a=9° case migrates 
over to and even crosses the centerline, but again the 
contours near the centerline are not characteristic of the 
presence of a vortex pair. 

Concluding Remarks 

Surface and flowfield measurements in a symmetric 
crossing shock/turbulent boundary layer interaction have 
been presented. From the results, the following conclu- 
sions can be made: 



z/<5 0 

Fig. 15 Volume fraction distributions for M = 
3.5, a = 2°, 6° and 9° cases, offset injection. 
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1) The length of the generator plate has a signifi- 
cant effect on the strength of the interaction for a given 
deflection angle. More specifically, the distance from 
the shock crossing location to the trailing edge expansion 
region must be large relative to the incoming boundary 
layer thickness in order to capture the true strength of 
the interaction. 

2) A qualitative difference is observed in the surface 
pressure distributions depending on whether or not a gross 
separation of the flow occurs. This feature may be useful 
in detecting flow separation when only surface pressure 
is available. 

3) Preliminary dynamic surface pressure measure- 
ments indicate, like previous studies, a significant un- 
steady characteristic in the interaction region. The dom- 
inant frequencies are generally centered around 1 kHz. 
For the strongest interaction considered, a high amplitude 
unsteadiness that occurs at ~150 Hz was detected. The 
undisturbed boundary layer power spectrum indicate the 
need for more extensive documentation of the acoustic 
properties of the NASA Lewis lxl SWT. 

4) Pitot pressure measurements illustrate the devel- 
opment of the boundary layer through the interaction re- 
gion and should be useful for CFD validation. 
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